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ABSTRACT 

We are developing time- and energy-resolved near-IR/optical/UV photon detectors based 
on sharp superconducting-normal transition edges in thin films. We report observations of the 
Crab pulsar made during prototype testing at the McDonald 2.7m telescope with a fiber-coupled 
transition-edge sensor (TES) system. These data show substantial (5a 0.3), rapid variations 
in the spectral index through the pulse profile, with a strong phase-varying IR break across our 
energy band. These variations correlate with X-ray spectral variations, but no single synchrotron 
population can account for the full Spectral Energy Distribution (SED). We also describe test 
spectrophotopolarimetry observations probing the energy dependence of the polarization sweep; 
this may provide a new key to understanding the radiating particle population. 

Subject headings: instrumentation: detectors, pulsars: individual PSR B0531+21 



1. Introduction 

Transition-Edge Sensor (TES) detectors are 
showing great promise as fast bolometer arrays in 
astronomical applications from the sub-millimeter 
through the X-rays (WoUman, et a/. 2000). In 
the near-IR/optical/UV range these devices offer 
good broad-band quantum efficiency (QE), high 
time resolution and modest energy resolution and 
saturation count-rate (Cabrera et al. 1998, Ro- 
mani et al. 1999). Along with competing cryogenic 
technologies (e.g. Superconducting Tunnel Junc- 
tion devices, STJs; Ferryman, Foden & Feacock 
1993, Ferryman et al. 1999) these sensor arrays 
offer the potential for important new capabilities, 
particularly for the study of rapidly varying com- 
pact object systems. We are developing sensitive 
TES spectrophotometer systems for such applica- 
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tions and report here on test observations of the 
Crab pulsar with a prototype fiber-coupled array. 

The basic principles and present performance of 
thin-film tungsten (W) TES devices have been re- 
cently summarized in Miller et al. (2000). In brief, 
the systems routinely achieve energy resolutions of 
^ O.lSeV at ^ 3eV, photon arrival time resolution 
of ~ 300ns and single pixel count rates of ^ 30kIIz 
without serious pile- up problems. En route to a 
high sensitivity, general purpose camera suitable 
for faint object high speed spectrophotometry, we 
have scheduled a number of astronomical demon- 
strations. Observations at the McDonald Observa- 
tory 2.7m Harlan J. Smith telescope were made in 
February 2000 with a system that incorporated a 
number of substantial improvements over the sin- 
gle pixel, fiber-coupled system used for the first 
astronomical observations described in Romani, 
et al. 1999. 



2. Experimental Apparatus 

The detector system used in these observations 
employs a 6 x 6 array of W TES pixels on a Si 
substrate with a 23/xm pitch, fabricated by our 
group using the Stanford Nanofabrication Facil- 
ity. For ease of fabrication l//rn wide Al leads, 
the 'voltage-bias rails', were connected in the pixel 
plane which means that the fill factor decreases 
towards the sides of the array as the active tung- 
sten area is reduced to accommodate these leads 
(Figure 1). The rails, superconducting during op- 
eration, are covered by an extension of the active 
W thin film of the TES pixel, providing electri- 
cal connectivity. Rails and pixels arc separated by 
Ifim gaps. Because of the underlying Al, photons 
absorbed above these rails in fact couple a larger 
fraction of the the deposited energy to the W e~ 
system than those absorbed directly on the adja- 
cent TES. These 'rail hits' thus increase the array 
fill fraction, but produce a 'satellite peak' in the 
energy PSF, which complicates the spectral anal- 
ysis (§3). The data system used for these tests 
had six read channels so in practice a 2 x 3 sub- 
array along the midline was used. These uniform 
20/xmx20/xm pixels gave an effective fill factor of 
96%. The geometrical peculiarities of the arrays 
used in these measurements are not intrinsic to 
the TES system. New arrays have masks to elim- 
inate the rail events. Buried wiring and focusing 
collimators have been designed that can further 
maintain uniform pixel size and high effective fill 
factors over substantially larger TES arrays. 

In operation, the Si substrate is cooled well 
below the ~ lOOmK W transition temperature, 
while the TES sensor is biased in a circuit with 
a fixed voltage. The resulting current allows the 
W electron system to self-heat into the middle of 
the superconducting-normal transition. The TES 
sensors are thus operating in a regime of strong 
negative electrothermal feedback (Irwin 1995). A 
major advance of the present apparatus is the first 
implementation of a NIST-dcvelopcd digital feed- 
back system which simultaneously linearizes the 
SQUID ammeters and functions as the data acqui- 
sition system. The digitized feedback signal mon- 
itors the TES device current. The data stream is 
processed through an FPGA with adjustable peak 
shaping and baseline-restore algorithms. In addi- 
tion to accurate GPS-stamped arrival times and 




Fig. 1. TES Array image. For these observa- 
tions, we recorded data from the 4 pixels along 
the midline best illuminated by the focussed fiber. 
The outermost pixels of the 8x8 configuration are 
not wired and serve only to mask the substrate. 

pulse height measurements for each photon, the 
system provides a number of quality control in- 
dicators including adjustable pile-up flags. Four 
of the pixels were read with independent digital 
channels of this new design; these provide the pho- 
ton set analyzed in this paper. This system is 
directly scalable to large numbers of read chan- 
nels and can be modified to allow multiplex read- 
ing of several pixels per read channel (Chervenak 
et al. 1999). 

There are two serious challenges in coupling 
these TES detectors to a typical ground-based as- 
tronomical telescope. One is to filter the incident 
beam so that the large flux of thermal photons (at 
A > 2/xm) from the warm optics does not satu- 
rate the TES sensor system, which has a modest 
maximum count rate. The second challenge is to 
maximally couple the beam of a large telescope 
to our present small (effectively 46/im) detector 
array. In the system described here we used a 
focal reducing train of flber and imaging optics 
(Figure 2). Given the relatively large plate scale 
(232/xm/arcsec at the f/17.7 Cassegrain focus of 
the 2.7m telescope) and best expected image of 
1-1.5" FWHM we adopted a 400/im (1.7 arcsec) 
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entrance aperture and fed the light through a 2:1 
reducing taper. This taper, along with the ~ 25m 
run of 200/xm-core fiber to the refrigerator bulk- 
head, were fabricated from low OH, low impurity 
(Polymicro FSU) Si/Si step-index fibers, providing 
both good IR and UV transmission. At the bulk- 
head of the cryostat, a portable "^Hc/^Hc dilution 
refrigerator, an ST-connector vacuum feedthrough 
passed the light to a ~ 5m length of high OH 
'wet' 200/im core fiber. The bulk of this fiber 
was spooled in thermal contact with the IK stage 
of the cryostat. In this way the strong OH ab- 
sorption bands of the wet fiber (Humbach 1996) 
acted to provide a cold filter with effective block- 
ing beyond 1.7/xm and several strong shorter wave- 
length absorption bands (e.g. at 0.95/im). These 
absorptions correspond to wavelengths of strong 
atmospheric OH emission. The system thus pro- 
vided transmission in the astronomically interest- 
ing J and H bands as well as reasonable through- 
put to 400nm where scattering in the fibers began 
to compromise transmission. 




Dilution refiigerator 
35 niKbase 



Fig. 2. — Schematic of TES Spectrophotometer 
system 

To illuminate the sensor array at the cold stage, 
we collimated the fiber beam with a 1.8mm 1/4 
pitch gradient index (GRIN) lens and then fo- 
cused with a 1.0mm diameter spherical lens to 
provide the best coupling to the central 2x2- 
pixel section of the TES array. Despite inevitable 
focal ratio degradation in the fiber train, the fi- 
nal //I focus with modest chromatic aberra- 
tion ensured that we had reasonable coupling into 



our central 160/im^ of active W TES. We mea- 
sured the absolute quantum efficiency of the sys- 
tem from the cryostat bulkhead to final recorded 
photons, including the ^ 50% absorptivity of the 
grey tungsten surface, using a calibrated 850nm 
light source, calibrated attenuators and an ab- 
solute power meter. These measurements deter- 
mined a total 850nm QE of 0.097. This was con- 
firmed by comparison with CCD integrations, us- 
ing cameras of known QE. The 'rail' hit photons 
contribute an additional 0.01 to the system QE. 
Inevitably, the fiber train to the telescope incurs 
additional losses. Unfortunately, these were ex- 
acerbated by a defect in the ST couplers of the 
taper section. As a consequence, the fiber section 
imposed a factor of « 2.5 decrease in the peak 
throughput, relative to the back end of the tele- 
scope. 

3. Calibration 

The energy scale for each pixel was set by inter- 
rupting the observations every few hours to take 
short ~ 60s exposures of a calibration lamp passed 
through a grating monochrometcr and fiber fed to 
the TES array. Typically we integrated with the 
1^* order set at ~0.5eV, l.OeV and 2.0eV. Observa- 
tions of the order peaks provided a quick calibra- 
tion of the system energy scale and non-linearity. 
In addition, to monitor any slow drifts in the en- 
ergy gain as might occur due to baseline temper- 
ature drifts in the dilution refrigerator, a series of 
calibrated 'heat pulses' with a range of energies 
were introduced to the W e~ system on a regu- 
lar cadence throughout the data exposures. These 
pulses produce photon-like signals of known ar- 
rival time and amplitude interleaved through the 
data stream. These heat pulses were excised and 
monitored for any sign of gain drift. In prac- 
tice, operation was sufficiently stable that no de- 
tectable drifts were generally seen in the heat pulse 
monitor. The monochrometer calibration spectra 
were used to form the response matrices Rk{Ei, cj) 
across detector energy bins Cj to the incoming pho- 
ton energies Ei for each pixel k. 

To complete the spectrophotometric calibra- 
tion, we observed the Bruzual-Persson-Gunn- 
Stryker spectrophotometric standard BD-l-30 2344 
at low airmass during each observing night of 
the run. Although convenient for its well cal- 
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Fig. 3. — The estimated total quantum efficiency 
of our system (fiber+detector) at the back end 
of the telescope. Atmospheric absorption, tele- 
scope inefficiencies and aperture losses give a typ- 
ical peak effective efficiency (relative to flux at the 
top of the atmosphere) of ^ 1%. For comparison 
the count spectrum observed from a 2.06eV set- 
ting of a grating monochrometer is shown. The 
low energy events (dashed line) include a rising 
continuum due to substrate (Si) hits, as well as 
a thermal IR H-band peak from the warm optics 
and lamp. 

ibrated absolute IR-optical-UV spectral energy 
distributions, at my = 10.2 the star had to be ob- 
served with substantial de-focus to avoid detector 
non-linearities at high count rates. We therefore 
also observed the my = 14.9 blue KPNO flux 
standard PG 0939+262 (Massey et aL 1988) to 
allow an improved estimate of the system UV 
response and a check on the point-source total 
efficiency, including aperture losses. The fiber fil- 
ter has strong absorption bands whose width is 
comparable to or smaller than our energy reso- 
lution. Also for steep red spectra, the satellite 
'rail hit' peak at 1.26 times the main peak en- 
ergy contributes significantly to the count rate. 
This made it impossible to extract the flux scale 
by directly dividing the predicted standard star 
counts by the observed channel counts. Instead 
the calibration star spectrum was convolved with 
the response function, the fiber absorption curve 



measured at high resolution, and the W absorp- 
tion spectrum to produce model calibrator count 
spectra. Comparing in the data space, an ef- 
fective area function varying slowly with energy 
was fit to model the additional focus and absorp- 
tion losses compared to stellar fluxes at the top 
of the atmosphere. Combining these factors, we 
derive model efficiency functions Ek{ej) for each 
active pixel k. After summing these we have a 
total (dctector+telescopc+atmosphcrc) efficiency 
for the system; this peaks at ~ 1.2%. In Figure 
3 we correct this efficiency for the estimated at- 
mospheric extinction, telescope inefficiencies and 
aperture losses to produce the effective response 
of the system (fiber optics -|- detector) at the back 
end of the telescope. For comparison, the en- 
ergy ranges corresponding to Johnson broad-band 
colors are indicated. Variation between the cali- 
bration sources indicate that the systematic errors 
in the flux scale are < 10% for the range 1-2. 5eV. 
In the blue, the calibration spectra have lower S/N 
and larger contribution from the rail hit events, so 
the systematic errors in the flux scale rise, prob- 
ably reaching 30% above 3.5eV. Comparison with 
archival Crab data below support the idea that 
our flux scale is systematically high by approxi- 
mately this amount in the U band. Imprecision in 
the estimates of the telescope and aperture losses 
also produce uncertainties exceeding 10% in the 
absolute flux scale. 

4. Crab Pulsar Spectrophotometry 

We observed a number of compact object 
sources, both accretion powered binaries and iso- 
lated spin-powcred pulsars. The relatively bright 
{rriy ^ 16.6) Crab pulsar was used to tune and 
characterize the system performance. 

Our GPS-stamped Crab photon arrival times 
were reduced to the solar system barycenter em- 
ploying the widely used TEMPO pulsar timing 
package (from pulsar.princeton.edu/tempo). We 
used the monthly Crab pulsar radio ephemeris 
published by Jodrell Bank (Lyne, Pritchard & 
Roberts 1999) to phase the photon arrival times. 
The light curves from our four nights of observa- 
tion are phased within ~ 10/is; the absolute phase 
agrees with the predicted infinite frequency pulse 
arrival time (referenced to the peak of the main 
pulse) to within ~ 30/U,s. This absolute precision 
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is rather better than the nominal ^ 100/is aecu- 
racy of the pubhshed monthly ephemeris and is 
completely adequate for our purposes, allowing a 
detailed comparison of our phase bin spectra with 
other (e.g. space-based) data sets. We focus here 
on a relatively long (3500s) coherent data set taken 
on February 6, 2000 (MJD 51581) for which the 
system parameters were stable and the observing 
efficiency was high. This data set provides much of 
our total S /N for the Crab pulsar and represents a 
unique time- and energy-resolved measurement of 
this well-studied object. A sample optical (1.5eV- 
3.5eV) light curve is shown for reference in Figure 
4. 
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Pulsar phase 

Fig. 4.— A sample 1.5eV-3.5eV (3550 - 8250 A) 
10^-bin light curve from 3500s of TES Crab data. 
Phase is referenced to the main peak of the ra- 
dio ephemeris. The upper histogram shows the 
same data at 100 bins (zero suppressed), magni- 
fied to show the baseline variation. Our selected 
background region is shown (label BG). 

There is appreciable structure in our light curve 
baseline even at pulse minimum. Previous stud- 
ies also infer magnetospheric emission throughout 
the pulse period (e.g. Sanwal et al. 1998). Re- 
cently Golden, Shearer & Beskin (2000) have ar- 
gued from MAMA imaging data that the phase 
interval 0.75-0.825 represents the DC pulse min- 
imum with ~ 0.4% of the integrated pulse emis- 
sion. Our data show significant structure in this 
region, but given the substantial background in 



our our 1.7" aperture, we select the phase inter- 
val (j)BG = 0.76 — 0.87 for our 'background' spec- 
trum, recognizing that there will be a small < 1% 
magnetospheric contamination; this will be domi- 
nated by Poisson fluctuations in nearly all of our 
extracted spectra. 
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Fig. 5. Phase- averaged TES Crab spectrum 
(full line with error flags) compared with broad- 
band colors. Predicted (extincted = 0.11; line 
histogram) and observed (filled dots) count spec- 
tra are shown for comparison (Right scale). 

In Figure 5 we show the phase-averaged spec- 
trum after background subtraction, convolution 
through the response matrix and application of the 
flux calibration solution of §3. To maximize the 
independence of the spectral bins, we first com- 
pute the expected count spectrum of the best-fit 
simple absorbed power-law model and then derive 
the Crab fluxes by multiplying the model flux by 
the ratio of the observed to predicted counts for 
each channel. This 'pre-whitening' avoids spread- 
ing fluctuations in the observed count bins be- 
tween final spectral bins and ensures that sen- 
sitivity to errors in the response matrices are of 
second order. Nonetheless, the derived spectrum 
represents only one model (the smoothest) that 
fits the count data. The error bars shown are 
Poisson only, derived from the independent counts 
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in the energy bins. For comparison, we plot op- 
tical/IR phase-averaged broad-band photometry 
tabulated in Eikenberry et al. (1997). We also 
show the phase-averaged spectrum fit by SoUer- 
man et al. (2000) to optical/UV data as the solid 
histogram - this is an = 0.11 powerlaw sub- 
ject to E(B-V)=0.52, R=3.1 interstellar redden- 
ing - displayed as a predicted count histogram. 
This may be compared with the observed spec- 
trum (solid dots). Notice that, in agreement with 
the IR photometry of Eikenberry et al. (1997), this 
powerlaw under-predicts the J/H flux. 

It is evident that our derived spectrum shows a 
somewhat steeper (redder) optical spectrum than 
the SoUerman et al. power-law fit and also lies be- 
low the B & U colors of Percival et al. (1993). 
This supports the suspicion of a systematic over- 
estimate of our sensitivity function above 2.5eV 
(possibly due to incomplete modeling of the rail- 
hit component of the PSF and/or pile-up). The 
notch in the spectrum at ~ 0.9eV is also likely 
an artefact caused by the rapid variation in the 
system sensitivity at the OH absorption band. 
Interestingly, we find that our E > leV spec- 
trum is however very well modeled by an extincted 
a = —0.2 power law, which is actually the same 
spectrum found originally by Oke (1969). The 
true spectral index is evidently quite sensitive to 
the reddening parameters and energy range chosen 
and, indeed, the optical/UV spectrum of SoUer- 
man et al. shows appreciable curvature. However, 
given the mismatch with the broad-band photom- 
etry in the blue, slight flux calibration systematics 
are the likely culprit. 

We can, of course, study spectral variations in 
the Crab by assuming that the phase-averaged 
spectrum is smooth and using our observed to- 
tal pulsar spectrum as the flux calibrator. This 
differential measurement allows a much improved 
removal of small scale variations in the sensitiv- 
ity function and provides robust detection of any 
spectral phase variations. Figure 6 shows a selec- 
tion of phase bin spectra, referenced to the total 
(phase-averaged) spectrum. There are evidently 
subtle but continuous variations in the broad-band 
spectral index and a distinct change in the be- 
havior in the near-IR band (see, for example the 
break at > leV in the leading edge of the main 
pulse, spectrum 1). The spectra show little fine 
structure except at the strong fiber absorption 



bands. Smooth continuum variations dominate, 
as expected for a broad-band emission process 
such as optical synchrotron. Errors are computed 
by propagating Poisson fluctuations in both the 
background subtracted spectra and the reference 
(phase-averaged) spectrum. 
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Fig. 6. Phase bin count spectra. The log of the 
photon flux is plotted against \og{E^) for several 
phase bins from the 13 contiguous bins indicated 
in the panel to the right. 

We can measure these smooth variations in 
smaller (p bins by fitting. Figure 7 shows a double 
power-law fit with a fixed break at 1.3eV. To ob- 
tain these spectral indices, we correct by a mean 
(phase-averaged) powerlaw. In view of the calibra- 
tion uncertainties above, we have chosen to refer- 
ence the spectra using a phase-averaged index of 
a = — 0.1. The spectral index fit errors are derived 
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from the Poisson spectral bin errors above. 
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Fig. 7. — Spectral index phase variations. Fits to 
optical (triangles, Upper panel) and near-IR (open 
circles, Lower panel) bands are shown. Note the 
rapid, anti-correlated variation in the main peak, 
indicating a varying spectral break. For com- 
parison, an optical light curve and X-ray RXTE 
(Pravdo, Angelini & Harding 1997) phase-resolved 
spectral indices are shown in the lower panel. 

The data show a rapid sweep in the spectral 
index, especially in the main peak. An optical 
light curve and peak phase marks are provided for 
reference. At the leading edge of the main peak 
the IR index does not reverse, but the spectrum 
becomes increasingly blue. This indicates strong 
break or depression of the IR at the leading edge 
of the main pulse. This is the origin of the IR 
variation in the leading pulse half-width noted by 
Eikenberry, et al. (1997). While the phase bins 
and error flags arc large in the bridge between 
the pulses, the optical and IR indices also differ 
strongly in this region, with a steep red IR spec- 
trum passing to a much harder component in the 
blue. This increase of the bridge flux at high en- 
ergies is in accord with the strengthening of the 
bridge emission in the X-ray and 7-ray bands. 

Phase-dependent spectral indices were first no- 



ticed in the Crab in the X-ray band (Pravdo & Ser- 
lemitsos 1981). The reference Rossi X-ray Tim- 
ing Explorer (RXTE) data above and more exten- 
sive BeppoSAX opscrvations reported by Massaro 
et al. (2001) show rapid variation in the hard X- 
ray power-law spectrum, especially near the first 
peak. The phase coincidence of the optical/IR 
variations imply a common origin, but there is 
evidently no simple relation between the indices 
in the two band. At higher 7-ray energies limited 
count rates do not allow the fine phase bins used in 
the X-ray analysis, but a correlated spectral vari- 
ation is seen across the Crab pulse profile (Fierro 
1995). 

In the optical, many authors starting with 
Muncaster & Cocke (1972) have used used time- 
resolved color photometry to suggest that the 
bridge region is bluer than the pulses. This is seen 
most clearly in the largely unpublished UBVR 
observations reported in Sanwal et al. (1998). Re- 
cently, Golden et al. (2000) have reported a red- 
dening in the bridge region, based on UBV pho- 
tometry, but this does not accord with other recent 
measurements. Note that with the exception of 
the Sanwal et al. (1998) data, limited count statis- 
tics forced all observations to use very coarse phase 
bins, typically summing over the entire pulse or in- 
terpulse. The a sweeps we see are locatcid largely 
within these components; this is a likely expla- 
nation for the non-detection of spectral variation 
in several studies (e.g. Carraminana, Cadez & 
Zwitter 2000, Ferryman et al. 1999). Thus the op- 
tical a variations through the pulses are seen here 
clearly for the first time. 

The simultaneous IR measurements of this data 
set also provide a new window on the pulsar 
physics, given the rapid variation in the spectral 
break below ~ IcV. Simultaneous optical/IR spec- 
tra were first reported in Romani et al. (1999). Re- 
cent simultaneous dual band (V & H) light curves 
have also been reported by Moon et al. (2001), 
but the low efficiency and limited S/N of these 
test data likely preclude measurement of the sub- 
tle spectral variations. Although some interest- 
ing indications of spectral variations can clearly 
be seen in the broad-band (non-simultaneous) IR 
photometry of Eikenberry et al. (1997), the data 
reported above are the first to show these fine- 
resolution spectral variations extending into the 
infrared. While the variations are highly signifi- 
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cant, longer integrations with better S/N arc cer- 
tainly desirable. Nonetheless, these data already 
provide some constraints on the radiation pro- 
cesses, which we describe briefly in §7. 

5. Crab Photon Statistics 

At our roughly kilohertz count rates, wo de- 
tected ~25 photons in each main pulse, giving 
us substantial sensitivity to photon correlations. 
In the radio most pulsars exhibit a Gaussian or 
log- normal distribution of pulse intensities. The 
Crab, however, also produces occasional radio 
pulses with up to ~ 10^ x the mean pulse en- 
ergy. These 'Giant pulses' are short duration 
bursts of flux, appearing randomly in narrow win- 
dows of pulsar phase. Lundgren et al. (1995) have 
searched for and placed limits on any correlation 
between giant radio pulses and 7-ray emission, 
while Patt et al. (1999) have found that no gi- 
ant pulses are present in the X-ray band. Percival 
et al. (1993) have also found no non-Poisson pulse- 
to-pulse variations in the optical flux of the indi- 
vidual pulse components. Lacking simultaneous 
radio data, we examined our data for any short- 
time scale photon correlations. 

We searched for non-Poisson photon distribu- 
tions by predicting the expected counts in a range 
of lag bins following each observed photon. Be- 
cause each detector has a small ~ lO/xs dead time 
associated with the recovery from a photon pulse, 
we searched the cross-correlations between our 
pixels, allowing us to extend the analysis to mi- 
crosecond times. To test the null Poisson hypoth- 
esis, we need an accurate estimate of the expected 
count rate in the face of the varying pulsed emis- 
sion and background in our aperture. This was 
obtained by forming a 100-bin light curve for each 
second of data and fitting the counts to a back- 
ground and pulsed flux amplitude. For each pho- 
ton arriving during this second, we used the two 
flt amplitudes to predict the photon arrival rate in 
the other three detectors on a range of timescales. 
The model flux in each lag bin was determined 
by an integration through the appropriate phases 
of a high-resolution (lO'^bin Figure 4) light curve 
from all the data in a selected energy range, scaled 
to the instantaneous detected flux. The observed 
counts can be compared to the Poisson expecta- 
tion for photons of a given energy range on vari- 



ous timescales and over various windows of pulsar 
phase. No strong energy dependency was noted, 
so we show the statistics of the full low background 
l-4eV energy range. 




PulsBr phase 

Fig. 8. — Ratio of average rate to Poisson rate for 
> IciV. The; grey scale image shows the ratio 
in independent lag bins, spaced as even octaves 
(vertically), vs. pulse phase (horizontally). Light 
shows high correlation and dark shows sub-Poisson 
counts. Bins with < 25 counts are set to 1. 

Figure 8 shows as a gray intensity scale the ob- 
served count rate in units of the Poisson predicted 
rate in 100 bins of pulsar phase. The vertical axis 
of the grey scale shows correlation bins covering 
time lags 6t = tsamp{2^ - 2^''^), where tgamp = 
0.31/is. runs from 1 to 11 (t = 0.3/iS to 635/is) 
in the gray scale plot shown. The lag bins are thus 
completely statistically independent. A signifi- 
cant correlation (or anti-correlation) on short time 
scales would be expected to show as a coherent 
bright (or dark) region at the appropriate phase, 
fading to grey in longer timescales. The data show 
fiuctuations in the individual amplitudes when the 
counts are low. There is a ~ 10 — 15% enhance- 
ment on timescales t ^ 5 — lO^s leading the first 
peak by ~ 1ms; this is of ~ 3a significance. There 
is also a ~ 30% enhancement at r < 3/is leading 
the second peak, but this is only ^ 2(t in the most 
significant lag bin. A few percent decrement at lag 
N = 7 appears correlated with noise in the 33/Lis 
bins in the fine template. All other significant fea- 
tures are smaller than 1%. We conclude that our 
photon arrival times are consistent with Poisson 
and that there is no evidence for strong correla- 
tions on short timescales, such as giant pulses, at 
optical energies. 

6. Crab Spectrophotopolarimetry 

Ultimately, complete exploitation of astronom- 
ical fluxes requires time-resolved spectrophotopo- 
larimetry. The strong, rapidly varying polariza- 
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tion of the Crab pulsar makes it an attractive 
test target. Phase resolved polarimetry has been 
presented in the optical (Smith et al. 1988) and 
the near UV (Smith et aL 1996), which give a 
basis for comparison. The double position-angle 
sweep reported in these papers has been shown 
to bo a good quantitative match to a rclativis- 
tic outer magnetosphere version of the rotating 
vector model by Romani and Yadigaroglu (1995). 
In this model emission from a range of altitudes 
forms caustics in photon arrival phase that define 
the pulses. Further, a combination of two distinct 
emitting regions can contribute to the intcrpulsc 
phases. If the radiating populations (and hence 
spectral indices and breaks) differ between these 
regions we would expect energy dependence of the 
polarization fraction. 

As a test, we observed the Crab pulsar on 
February 7, 2000 (MJD 51582) with a simple po- 
laroid filter placed before the focal plane. We 
integrated at two full cycles of position angle 
with the following exposures: 0° (1370s,630s), 45° 
(520s,690s), 90° (740s,620s), 135° (960s,680s) with 
angles referenced to the equatorial coordinate sys- 
tem. Of course, varying background and vary- 
ing aperture losses between the exposures mean 
that the light curves must be renormalized be- 
fore forming q = {Iq - hoj/ih + -^90 ) and u = 
{hb — Ii35)/{l45 + -^135)- To subtract the sub- 
stantially polarized nebular emission in our aper- 
ture, wc employed our usual background phase 
window (0.76-0.87), following the prescription in 
(Smith et al. 1996). While residual pulsar emission 
at these phases may be highly polarized (Smith 
et aZ. 1988), the amplitude is small enough that 
its contamination of the peaks and bridge region 
should be smaller than our statistical errors. An 
examination of the Smith et al. (1988) data (kindly 
supplied in digital form by F.G. Smith) showed a 
polarization minimum during the bright first pulse 
in the region </> =0.005-0.025, consistent with 0, 
within the error bars. This provides a convenient 
exposure monitor based on the pulsed omission it- 
self. These data were taken in a broad optical 
band with the 'RGO people's photometer'; we syn- 
thesized this band in our data, selecting the appro- 
priate energy bins and using the polarization min- 
imum flux as a monitor of the total Crab exposure 
at each polarization position angle. 

In Figure 9, we show the resulting polarization 
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Fig. 9. — Test Crab Spectropolarimetry in three 
optical bands. Two pulse periods are shown for 
clarity. Above: Position angle. An optical light 
curve is included for phase reference. Below: Lin- 
ear polarization fraction. 

fraction p = (q^ + u'^Y^'^ corrected for the wave- 
length dependent polarization efficiency of the fil- 
ter (measured using monochrometer light passed 
through crossed filters). We also corrected for 
the bias due to non-normal statistics at small po- 
larization, with the observed polarization po cor- 
rected to p = po[l — (f/Po)^]^^^- In the up- 
per panel we show the polarization position an- 
gle = 1/2 x Arctan(u/g). The error flags shown 
are the propagated Poisson errors. Differences be- 
tween the two sets of position angle data were 
within the statistical errors, so these were summed 
in the polarization amplitude and sweep shown. 
Limited statistics restricted us to 18 phase bins 
through the pulse, so fine structure is likely not 
fully resolved. 

The overall polarization behavior follows that 
measured in previous optical/UV experiments 
quite closely. We show the optical data divided 
into three energy bands. The rapid decrease in 
polarization and strong position angle sweep of 
the first pulse appear quite energy independent. 
There are small differences in the polarization and 
position angle of the bridge and second pulse, but 
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these are of marginal significance. Given that the 
phase variations in the spectral index and break 
are strongest below 1.5eV (§4) we might expect 
the strongest spectral effects in the near-IR. In- 
deed, our near IR data do appear to show signif- 
icantly different polarization properties. Unfortu- 
nately, the increased IR background and difficulty 
in calibrating the rapidly varying IR polarization 
response of the filter preclude a quantitative com- 
parison with the optical data in the present exper- 
iment. Serious phase-resolved spectropolarimetry 
will require imaging solutions for local background 
subtraction; as part of our next-generation imag- 
ing system (see below), we plan to install Wollas- 
ton prisms for imaging spectrophotopolarimetry 
over reduced fields. 

7. Spectral Variations and Conclusions 

The phase variations through the Crab pro- 
file described above clearly provide new tests of 
the physics of the optical/IR radiation. The po- 
larized IR through hard X-ray non-thermal emis- 
sion from the Crab is generally held to be syn- 
chrotron radiation from a population of prop- 
agating and cooling in the outer magnetosphere 
(e.g. Lyne and Smith 1990 and references therein; 
Romani 1996; Wang, Ruderman & Halpern 1998). 
The broad-band phase-averaged SED is smooth, 
with a varying from —0.7 (hard X-rays) through 
—0.5 (soft X-rays) to ~ at optical frequencies. 
There is evidently a peak in the near-IR and then 
the spectrum drops off dramatically at lower en- 
ergies (Middleditch, Pennypacker & Burns 1983), 
although this important result definitely requires 
confirmation. 

Recently, Crusius-Waetzel et al. (2001) have 
provided a synchrotron interpretation that pro- 
vides a plausible explanation of the soft-X-ray 
through IR phase averaged SED in terms of a 
single N oc E'^^ population (presumably pro- 
duced through an pair cascade in the outer 
magnetosphere) with s = 2. In the X-ray band, 
the emission is assumed to be classical relativistic 
synchrotron radiation with Fe oc E~^^~^^/'^; this 
gives the observed a ~ —0.5 spectrum observed 
for the powerlaw X-ray emission of the Crab and 
other young pulsars. However, Crusius-Waetzel 
et al. additionally argue that the population 
should have a component with very low pitch an- 



gle ^ <^ 1/7, for which the synchrotron emission 
formulae give Fe oc E'^~^, when averaged over 
the radiating beam. For the Crab, they infer 
7 ^ 10^ - 10^ and * < 10^^, and connect the 
a « spectrum with that seen in the optical. Fi- 
nally, they suggest that below leV the particles 
dominating the emission spectrum have even lower 
7, so that the observed pulse width ~ 0.01 
samples only a subset of the I/7 radiation cone 
and argue that this leads to a Fe cx iJ"*"* spec- 
trum. This index is taken to be, at least in part, 
an explanation of the steep rise in the mid-IR. 
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Fig. 10. — Low pitch-angle prediction for the op- 
tical spectral index variation (full line) compared 
with optical and IR TES data. 

If wc take the X-ray spectral index to be a lo- 
cal measurement of the electron powerlaw s then 
this picture predicts the correlated optical spec- 
trum in the low pitch angle limit. In Figure 10, 
we show the TES data compared with this pre- 
diction, based on the 0. 1-4.0 keV data of Mas- 
saro et aZ. (2001). Clearly this phase bin appli- 
cation of the low pitch-angle picture is not ade- 
quate: while the index and trend for the second 
pulse are reasonably well predicted, the spectra of 
the bridge region and, especially, the first pulse 
show no match to the extrapolations of soft X- 
ray data. When one considers that the spectral 
index in all phase bins steepen into the hard X- 
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rays, then one must conclude that a simple ap- 
plication of powerlaw synchrotron emission, even 
with a phase- varying s, is not adequate. 

One factor likely important in a more com- 
plete model is the inclusion of the local beam- 
ing of the radiation cones along the emission sur- 
face. Given the sharp pulses and the observation 
at a single line-of-sight, the radiation must differ 
from the emission-cone averaged synchrotron ex- 
pressions above. Models including these effects are 
sensitive to the detailed geometry. As an exam- 
ple, consider the strong IR break at the leading 
edge of the main pulse. In the outer magneto- 
sphere caustic picture of the pulse formation, the 
leading edge of the main pulse arises from emis- 
sion propagating along the radiating surface for 
the longest distance. Hence, one would expect 
this to be the region first showing self-absorption 
suppression. While the overall IR suppression of 
the main peak reported by Penny (1982) sugges- 
tive of uniform self- absorption is not seen in our 
data, the observed leading edge suppression and 
steepening may be a more restricted instance of 
self-absorption. Polarization measurements ex- 
tended into the IR will be a good discriminant 
between absorption effects and emission breaks, 
as a restricted sampling of the synchrotron emis- 
sion cone leads to polarization signatures varying 
significantly in accord with the spectral breaks. 

We are developing a focal plane camera based 
on TES detector arrays, with the potential for 
time resolved imaging spectral polarimetry. As 
the data of Golden et al. (2000) and Ferryman 
et al. (1999) show, the absolute photometry and 
backgroTind subtraction of such true imaging mea- 
surements can provide much improved isolation 
of faint pulsar phenomena. With the clean mo- 
noenergetic energy PSF of a 'masked' TES array, 
improved system QE and improved energy resolu- 
tion, observations of the Crab (and other fainter 
pulsars) should allow a detailed examination of 
the phase-resolved phenomena discovered in the 
exploratory measurements reported here. 
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